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. We discuss some issues related to the computation of black hole entropy in loop 

I quantum gravity from the novel point of view provided by the recent number- 

ed I theoretical methods introduced by the authors and their collaborators. In particular 
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we give exact expressions, in the form of integral transforms, for the black hole en- 



' tropy in terms of the area. We do this by following several approaches based both 

on our combinatorial techniques and also on functional equations similar to those 
employed by Meissner in his pioneering work on this subject. To put our results in 
O I perspective we compare them with those of Meissner. We will show how our methods 

confirm some of his findings, extend the validity of others, and correct some mis- 
^jj^' takes. At the end of the paper we will discuss the delicate issue of the asymptotics 

of black hole entropy. 
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I. INTRODUCTION 



The derivation of the Hawking-Bekenstein area law for realistic black holes is one of the 
most important achievements of loop quantum gravity (LQG) [l|, Isj. This, together with 
. the recent insights on the Big Bang singularity provided by loop quantum cosmology 
r> I are two of the physical pillars on which the formalism is currently supported. The fact 
I that areas are quantized in LQG is another result that plays a relevant role in the study of 
black hole physics. This is so because, as pointed out in classic papers by Bekenstein and 
Mukhanov |5| , there are good reasons to believe that black hole areas should be quantized in 
such a way that the spacing between consecutive values of them is constant. Here, however, 
LQG does not seem to provide a completely satisfactory description because the spectrum 
of the operator is not equally spaced. 

A surprising development in this matter took place when Corichi, Diaz-Polo, and 
Fernandez-Borja found, by using a direct computer intensive approach, that effectively 
the black hole degeneracy spectrum for small black holes can be considered to be equally 
spaced. This intriguing result -obtained for black holes of around a hundred Planck areas- 
is beautiful because it suggests that, after all, the expected behaviour for the entropy can be 
somehow obtained within the LQG framework. However appealing this finding may seem. 
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it should be taken with some care because the description used for black holes in LQG, 
modeled by isolated horizons (IH), can only be approximate^, though arguably good for 
sufficiently large objects. It is then very important to find out if the observed microscopic 
behaviour of the entropy is also present for macroscopic objects. This would be a very nice 
result because the Beckenstein-Mukhanov prediction would be non-trivially realized within 
LQG in a macroscopic regime where the model used to describe black holes is arguably ac- 
curate. The main obstacle to find out if this is the case is the impossibility of extending the 
numerical methods used to date to these large scales. Instead, one must rely on a suitable 
asymptotic analysis based on closed and explicit expressions for the black hole entropy. 

This problem has been already considered in the literature; in fact a solution to it is 
proposed in a well known paper by K. Meissner [3|. According to the result presented there 
the asymptotic behaviour of the entropy does not display the oscillations described in 
What is happening then? Is the behaviour observed in [q] an artifact of the algorithm used 
to compute the black hole entropy? Is there an independent way to check it? 

An answer for these questions appears in j3] where a new method to compute black hole 
degeneracies based on a number-theoretical approach is developed. The procedure proposed 
in that paper provides an algorithm that can be used to check and extend previous numerical 
results. The conclusion is unambiguous: the results on the entropy originally found in 
are quantitatively correct and persist up to areas an order of magnitude larger than the ones 
considered in that paper^. The structure in the entropy spectrum is clearly present in the 
new regime explored in ^] and the constancy of the effective spacing between the areas is 
confirmed. 

Once the reality of this effect has been settled beyond doubt it is natural to pose an 
additional set of questions, for example: is it possible to find a single procedure that allows 
us to derive both the microscopic results of 0, 0] and the expressions given by Meissner 
for the macroscopic black hole entropy? Can we trust the asymptotic analysis performed 
by Meissner that seems to exclude, in the large area regime, the behaviour found for small 
black holes? 

Some steps towards answering the first of these questions have been taken in where the 
new number-theoretical methods introduced in [7|] were used to obtain generating functions 
for the black hole degeneracy spectrum. This is an important starting point because once 
closed expressions for them are available it is possible to obtain exact expressions for the 
black hole entropy in terms of the area. The first goal of this paper is to explain how this 
can be done. 

After the exact expression of the black hole entropy is obtained we will compare it with 
the ones given in j^. As we will see, when the projection constraint is not taken into account, 
we exactly reproduce^ the result given by formula (13) of [3|. On the other hand when the 
projection constraint is incorporated our results disagree with those of Meissner. In order 
to understand the source of this discrepancy we solve the problem from scratch by using 
functional equations as done in [3|. This is the second goal of the paper, namely, rederive the 
results of Meissner by using the right functional equations. As we show there is an error in 



^ Notice, for example, that microscopic black holes should evaporate very quickly; a fact that cannot be 

easily taken into account with the techniques currently available in LQG. 
^ Actually up to the largest areas that we have been able to reach by using laptops or personal computers. 
^ In order to be fully consistent with the definition of entropy that we are using here a term should be 

added to formula (13) of [3]. We will comment on this later. 
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the relation used by him when the projection constraint is included. We identify the source 
of this error and correct it. Furthermore, the inclusion of a superfluous pref actor in the 
equation given in Q makes it difficult to solve in closed form. As we will see it is possible to 
write a simpler -but equivalent- functional equation and solve it exactly in a straightforward 
way. After doing this we get the same expression for the black hole entropy that we find 
by using our generating functions. In our opinion this provides a solid confirmation of our 
results and highlights the power of the combinatorial and number-theoretical methods of 
0, Isl. Along the way we show that, at variance with the claims of Meissner^, it is in fact 
possible to give exact expressions for the entropy, without using any approximation, also 
when the projection constraint is taken into account. 

After finding exact expressions for the black hole entropy as suitable integral transforms, 
we want to raise some points concerning their asymptotic behaviour and then answer the 
second question posed before. We do this only in the case where the projection constraint 
is not used (the difficult problem of obtaining the asymptotic behaviour of the entropy in 
full generality will be addressed elsewhere). Notice, however, that as shown in 0, [1], the 
interesting structure of the black hole entropy is present even if the projection constraint is 
not taken into account so it makes sense to concentrate on this somewhat simpler situation. 
The most important issue to discuss now has to do with the poles of the Laplace transform of 
the entropy. The main result that we prove in the paper in this regard is that the real parts 
of these poles have an accumulation point precisely for the value of the Immirzi parameter 
7m computed by Meissner. In order to show this we need some auxiliary results concerning 
the distribution of poles that we will explicitly write down and prove. The fact that the 
real parts accumulate to 7m makes the asymptotic study of the entropy highly non-trivial 
as we will discuss. The bottom line is that the behaviour of the entropy in terms of the 
horizon area A is somehow proportional to exp(7M^), however it is not possible to exclude 
the possibility that this behaviour is modulated by an oscillatory, possibly decaying, term. 
Understanding the details of this behaviour is a crucial issue in order to see if the structure 
found in the entropy spectrum for small black holes is present in the macroscopic regime. 
We want to stress at this point that, according to our results, there are good reasons to 
believe that the value of the Immirzi parameter jm computed by Meissner is actually the 
correct one. At any rate, a final statement on this fact can only be made when the full 
asymptotic behaviour of the entropy is found. 

The paper is organized as follows: After this introduction we start with section [TTl where 
we discuss in detail the obtention of black hole entropy in LQG. The presentation that we 
give here is complementary to the one appearing in where we gave a unified treatment for 
the different types of countings proposed in the literature. Here we will use only the standard 
entropy definition of black hole entropy in LQG as spelled out in We show in section 
mil how one can obtain exact formulas for the black hole entropy by using the generating 
functions appearing in Section [IV] considers the same problem by using functional 
equations in the spirit of Meissner. As we will show we recover the results obtained in 
section [ITl without resorting to any approximation. These are then compared to the ones 
obtained by Meissner in [3|. Section |V] is devoted to a study of the inverse Laplace transform 
giving the black hole entropy in the simplified case where no projection constraint is used. In 
particular we enunciate and prove several lemmas concerning the distribution of the poles 



^ According to Meissner his formula (35) is obtained after neglecting some contributions in the resolution 
of the functional equation that gives the entropy. 
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of the integrand. We will prove an important result on this: the real parts of the poles 
have an accumulation point precisely at 7m- This may invalidate the conclusion that the 
asymptotics of the entropy as a function of the area is proportional to expi^jM^) as will 
be shown with a concrete example. We end the paper in section IVII with our conclusions, 
comments, and a short review of the problems that remain to be solved in order to fully 
understand the macroscopic behaviour of black hole entropy within LQG. 

II. COMPUTATION OF THE ENTROPY IN LQG: THE BLACK HOLE 

DEGENERACY SPECTRUM 

This section describes the number-theoretical and combinatorial approaches of 0, @] to 
compute black hole entropy in LQG. Here we will consider only the standard counting of 
[2|. Quoting almost verbatim from this paper we take the following definition: 

Definition II. 1. The entropy S{a) of a quantum horizon of the classical area a, according 
to Quantum Geometry and the Ashtekar-Baez-Corichi-Krasnov framework fi], is 

S{a) = logn(a) , 

where n(a) is 1 plus the number of all the finite, arbitrarily long, sequences fh = (mi, . . . , m„) 
of non-zero half integers, such that the following equality and inequality are satisfied: 

n n 

J2mi = 0, J] V|m,|(|mi| + 1) < 

i=l i=l P 

Here 7 is the Immirzi parameter of Quantum Geometry and ip the Planck length. The extra 
term 1 above comes from the trivial sequence. 



Let us start by introducing some notations, unit conventions, and definitions. In the 
following 

N = {1,2,3,...}, No = {0,l,2,3,...}, Z = {0, ±1, ±2, • • ■ } , = {±1, ±2, ■ • • } . 
We will also define 

Z/2 := {0,±l/2,±l,±3/2,---} 

with analogous definitions for N/2, No/2, and Z^,/2. The Kronecker symbol is written as 
and 9{x) := X[o,oo)(^) denotes the step function that we use throughout the paper 
(the characteristic function of [0, cxd) satisfying 6'(0) = 1). Finally [a;J denotes the integer 
part (floor) of the real number x. In our previous work on this subject 0, 0] we have 
used units such that 47r7£p = 1. Here, however, we will take 87r7£p = 1 to facilitate the 
direct comparison of our results with those of Meissner. A simple translation guide between 
expressions written in the two different unit systems can be given: in order to take formulas 
from the 87r7£p = 1 unit system to the 47r7£p = 1 one, it is enough to substitute the areas 
a appearing in them by a/2. 
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It is convenient now to define several sets that will play a relevant role in the following. 
First, given {a,p) G [0, oo) x Z/2, let J\f<{a,p) and Af<{a) be the sets 

n n 

X<{a,p) := {m\3n eN : m e (Z*/2)", ^ ^/\mi\{\mi\ + 1) < a, = p} , 

i=l i=l 

n 

Af<{a) := {m\3nen : rh e (Z,/2)" , ^ ^/\mi\{\mi\ + 1) < a} , 

1=1 

and lets us denote by N<{a,p) and N<{a) their respective cardinalities. Notice that the 
entropy S{a) is given by^ 

e^(") = n(a) = A^<(a,0) + 1. (2.1) 
The set A/<(a) can be written as the disjoint union 

Ar<(a)= U Ar<(a,p) 

and hence the cardinality N<{a) can be obtained in terms of the numbers N<{a,p) by 
summing in p, i.e. 

iV<(a)= J] iV<(a,p). 

pez/2 

We will also consider the sets 

n n 

M{a,p) := {m\3nen : me (Z,/2)" , ^ V|mi|(|mi| + 1) = a, =p} , 

1=1 i=l 

n 

Ar(a) := {m I 3n G N : m G (Z,/2)'^ , ^ ^\mi\{\m,\ + 1) = a} = |J Ar(a,p) , 

whose cardinalities N{a,p) := |A/'(a,p)| and A^(a) := |7\/'(a)| satisfy 

iV(a)= 5^ N{a,p). 

pez/2 

It is clear that the AQ-sets can be written as disjoint unions of A/'-sets. Explicitly 
J\f<{a,p) = U Af{a',p) , Af<{a) = [j U{a') = IJ U ' 

a'<.a a'<a a'<apgZ/2 

and hence 

N^{a,p) = Y,N{a',p), iV<(a) = iV(a') = ^K'^)- 

a'<a a'<a a'<apgZ/2 



^ The number 1 introduced in this definition was not taken into account in [3| so the formulas appearing in 
that paper should be corrected accordingly. 
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Notice that in order to compute the black hole entropy according to the definition given 
above we only need to know iV<(a, 0). However, it is convenient at times to work with the 
N<{a,p) so we will keep the p-label in the following and impose the condition p = only 
when needed. 

We will next obtain exact formulas for N{a), N{a,p), N<{a), and N^{a,p). The first two, 
N{a) and N{a,p), refer to what we call the black hole degeneracy spectrum jl] whereas the 
last two are directly related to the black hole entropy. Before giving a formal derivation we 
will summarize the procedure that we will use. The first step in all the cases is determining 
the finite sequences of arbitrary length n consisting of non-zero, positive, half integers |mj|, 
i = 1, . . . ,n, satisfying the condition 

n 

^/\mi\{\mi\ + 1) = a 

i=l 

for a given value of a. This can be done by first finding out the possible values for the \mi\ 
(with their multiplicities) compatible with a and then considering all the distinguishable 
permutations of them. This first step basically solves the problem of obtaining N{a). To 
obtain A^(a,p) one must take into account the projection constraint 'rij = p. As 

explained in [9[] this can be easily done by using generating functions in the form of Laurent 
polynomials. Finally, to find N<{a) and N<{a,p) one has to give a method to add the 
cardinalities given by N{a') and N{a',p) for all the eigenvalues a' in the area spectrum 
smaller or equal to the given a. 

Let us start by considering the sets A/'(a) and M{a,p). We want to give generating 
functions for both N{a,p) and N{a). To this end, let us consider first the auxiliary set 

n 

K{a) := {A?|3n G N : G N" , ^^{ki + lf -I = 2a } . 

i=l 

Now, there exists a surjective map vr : Mi^a) /C(a) defined by 

A/'(a) n (Z*/2)" 3 m t-^ 7r(m) = A; G fC{a) n , where h := 2\mi\ , 2 = 1, . . . , n. 
Clearly, if m G A/'(a) and k = 7r(m) we get 

J2 V|m,|(|m,| + 1) = a^Yl Viki + l^-l = 2a . 

i i 

The map tt is not injective and hence given k G /C(a) fl N" it is not possible to unambigu- 
ously reconstruct rh because there are several m G Af{a) such that '7r(m) = k (i.e. several 
acceptable choices for the signs of the m,), 

rh G TT~^{k) <^mi e {-ki/2, ki/2} , i = l,...,n. 

In order to determine the cardinality K{a) = |/C(a)| let us fix a G [0, cxd) and consider the 
equation 

J2 Viki + iy-l = J2 + 1)2 - 1 = 2a , (2.2) 

i=l k=l 
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where denotes the number of times that the integer k E N appears in the sequence k. 
Equation ( 12. 2p should be understood as an equation in the set of unknowns {(/c, n^)}. It is 
important to reahze that once the set of all possible solutions is determined the sequences 
k can be found by considering all the permutations of a multiset where each k appears 
times. 

Notice that we can always write ■\/(A; + 1)^ — 1 as the product of an integer times the 
square root of a square-free positive integer number (SRSFN) by using its prime factor 
decomposition. Hence, K{a) > implies that a is constrained to have the form 



2a = J2^^VP^ (2-3) 



i=l 



where qi G No and Pi > I are square-free integers (we will consider that pi = 2 < p2 = 3 < 
■ ■ ■ < Pi < pi+i < ■ ■ ■ , etc) . In order to proceed from here we must first identify the allowed 
values of k such that ^J{k + 1Y — 1 is an integer multiple of some y/pl, and then determine 
the value of that tells us how many times each of them appears. We deal with the first 
problem by solving the Pell equations associated to each of the SRSFN's in the r.h.s. of 
O, i.e. 

^/{k + lf-l = y^, ^{k + If - p,y^ = 1, (2.4) 
with y E N. The solutions can be labeled as {{kl^,yU : a G N}, where here the index i 



refers to the square-free numbers in each of the Pell equations (see, for instance, jlOl for 
details on the Pell equation). Once these numbers are known the can be found [7| by 
solving the system of r-uncoupled, linear, diophantine equations 



^y'oJ^ki,=qh « = l,...,r. (2.5) 
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Notice that, once the qi are fixed, only a finite number of labels fcj,, a = 1, . . . , Mi{qi) := M, 
come into play in the equations (12.51) . It may happen that some of these equations admit 
no solutions. In this case 2a = Yll=i Qi\/Pi is such that K{a) = (i.e. a does not belong to 
the area spectrum). On the other hand, if they do admit solutions, 2a = ^^^^ Qiy/pl is such 
that K{a) > 0, the numbers fc^ tell us the allowed values of the components of k involved, 
and the n^i^ count the number of components of k whose values are equal to fc^. 
Let us now define the set 

r Mi 



Given k = (ki, ■ ■ ■ , k.,.) G K{a) it is possible to construct 
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sequences in /C(a) by permuting the elements of the multiset^ Ujkj. Then, given 2a = 
Yll=i liy/Vii the cardinahty of /C(a) is 

A (a) = |/C(a)| = — j- • 

llj=l ll(fc,nfe)eki "'s- 

In order to determine N{a) from K{a) it is enough to reahze that given k G /C(a) each 
component /cj gives rise to two different values rrii G {—ki/2, ki/2} of the corresponding 
components of m G 7r~^{k), whereas the number of allowed configurations once the projection 
constraint is incorporated can be easily obtained by using a simple generating function as 
in [9|. We have then the following 



Theorem II. 1. The value of N{a) is given by 

whereas, when the components of the sequences m are required to satisfy the projection 
constraint rrii = p, we get 

Nia,p):=[z^^] J2 ( n-"n "''' n n + ■ (2.7) 

□ 

We have used the notation [z'^^]f{z) for the coefficient that multiplies 2;^^ in the Laurent 
expansion of the function /. The numbers N{a,p) can be conveniently encoded in the 
generating function 

G{z;x„X2,...) = (l-5^5^(z'="+^-'=^)xf^)" . (2.8) 

i=l 0=1 

The coefficient of the term z'^^xf ■ ■ ■xf' tells us the value of N{a,p) for 2a = ql^/Pl + ■ — I- 
Qry/Pr- Moreover, using the relation N{a) = J2pez/2-^i^^P)y generating functions for 
the numbers N{a) can be obtained from the generating function for the N{a,p) numbers 
simply by setting the z-argument equal to one. Hence, the coefficient of the term xf^ ■ ■ ■ x^^ 
in the power expansion of 



G(l;xi,X2,...) = (l-25^5^xf^)"' (2.9) 



1=1 a=l 



^ Here we are considering the sets as a multisets. Hence (fc, rik) E should be understood as a notation 
signaling the presence of elements, each of them equal to fc, in the multiset. 
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gives us the value of N{a) for 2a = qi^/pi + ■ ■ ■ + (ir^fPr- It is important to notice that 
G(^z] 0, O^. . . ) = 1. This can be interpreted as the extra one that appears in the prescription 
fJ2.ip of [2[]. Formulas (12. 7p and fl2.6p summarize in a compact way the algorithms used in 
[3] to compute the black hole degeneracy spectrum. The same information is stored in the 
generating functions (12. 8p and (12. 9p in a way that will let us write down very compact exact 
formulas for N^{a,p) and N<{a). This is the purpose of the next section. 

III. FROM N{a),N{a, p) TO N< (a) , N< (a, p) 

The coefficients of the power terms of the generating functions given above can be writ- 
ten in closed form as multiple contour integrals depending on gi, . . . , by using Cauchy's 
theorem. For example, given 2a = qi^/pi + • ■ ■ + (ir^fPr , 

N{a) = [xf---x'^r]G{l;xi,X2,...,Xr,0,...) 

with appropriately chosen contours 7^ surrounding the origin. Such integral representations 
are usually a good starting point to obtain asymptotic expansions in terms of the parameters 
appearing in them (the qi in this case). As we are really interested in the asymptotic 
behaviour of the entropy as a function of the area there is a necessary intermediate step: 
determining the parameters qi{a) as functions of the area. If these can be written in a 
reasonably simple closed form and are regular enough, just by plugging them in integral 
expressions of the type written above we could have closed formulas in terms of the area 
for the objects that we are interested in. The unfortunate fact is that the coefficients qi{a) 
as functions of the area oscillate wildly and in an rather unpredictable way so this direct 
approach is, to say the least, hard. 

The sector = {a„ : n G N} of the spectrum of the area operator relevant in the study 
of black hole entropy is a countable, ordered {an < ctn+i for all n E N), subset of the real 
line. It is then possible, in principle, to build the sequence {A^(a„,p) : n G N}. For a fixed 
value of the area spectrum a„ we can then obtain N<{an,p) as 

n 

N<{an,p) = '^N{ai,p). 

i=l 

In fact, if the values of N^a^p) are encoded in the generating function gp{x) = 
^^gj^ A^(a„,p)x" this summation can be carried out by a well-known procedure consist- 
ing in taking the new generating function^ 

G,{x) = ^ = TN^{an,p)x-. 
1 — X ^-^ ~ 

Again this is difficult in the present case because to succeed with this approach one would 
need to have an appropriate way (i.e. a manageable closed formula) to find the numbers q^ 



^ For simplified models where the area spectrum is taken to be equally spaced this can be done 
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corresponding to the n*'^ element of the ordered set of area eigenvalues Am- This can be 
seen to be equivalent to solving the following two problems: 

i) Given an eigenvalue of the area a e Am, how many smaller eigenvalues do exist? (we 
refer to this as the area ordering problem). 

ii) Given ^4 e M, what are the values of the qi corresponding to the largest area eigenvalue 
a e satisfying a < A7 (alternatively to the closest eigenvalue to ^4). 

The last question, in particular, must be answered because in practice we want to treat the 
areas as a continuous real parameter. Although the problems just described are probably 
not unsurmountable, here we will avoid them and give a remarkably simple procedure to 
perform the required summations by using Laplace transforms^. The key idea is to make 
use of the following two simple facts: 

1) jC[5{a — a);s\ — e'"'^ for a > 0. Here 6{a — a) denotes the Dirac delta distribution 
supported at a. 

2) If C[F{a); s] = f{s), then the Laplace transform of F{a') da' is simply s~^f{s). 

Generating functions and Laplace transform: Let us consider the "generating func- 
tion" ^ 

G{x)=Y,Pnx"- (3.1) 

neN 

defined in terms of the sequences A — {an ■ n e N}, < < q;„+i, and B — : n e N}. 
Notice that, in general, the are not necessarily integers. If we take the distribution 

^(") = X]/5n5(a - an) 

neN 

we have that 




represents the sum of the numbers /3„ corresponding the the values of q;„ smaller than a 
(for a ^ A; notice that, otherwise, the integral in the above formula for F< is ill-defined). 
In order to extend the definition of F< to the values of a G ^ as the sum V] f„.^ Bn it 
suffices to consider the limit from the right 

F<{a) = lim / F{a')da' = V /3„ , a > 0. 

^-^«+ Jo I ^ X 

{n:an<a} 



In the following jC[F{a);s] denotes the Laplace transform, expressed in the variable s, of the function 
F{a). On the other hand, jC~^[f{s); a] denotes the inverse Laplace transform of the function f{s) in terms 
of the variable a. 

We will not worry about regularity or convergence issues here. 
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If /3n ^ for every n G N the existence of these hmits is guaranteed by the fact that F< is an 
increasing function. If not all the values of are equal this function has jump singularities 
in the values a„. It is important to realize that the spacing between the values of a„ plays 
no role in the previous formula, so it works equally well for evenly or unevenly spaced values 
of a„. Also it is important to realize that under mild conditions on the sequences A and B 
the function '^{n-an<a} Pn0{ci — C(n) will be exponentially bounded and its Laplace transform 
well defined. 

In order to compute F< (a) the idea is then to consider the function 

P(s) = G(e-^) = 5^/?„e-°''^ 

neN 

defined in terms of (13.11) . Notice that it encodes the information about both sequences A 
and B. We can think about P{s) as the Laplace transform 

C[^Pn5{a - a„); s] = P{s). 

neN 

The arguments given above immediately tell us that 

F^ia) = C-'[s-' Pis); a] = C-'[s-'J2f3ne-''-' -,01 

neN 

if a does not coincide with any of the values corresponding to the sequence A. On the 
other hand if a G ^ then the fact that at jump singularities the inverse Laplace transform 
always gives the average between the left and right limits means that C~^[s~^P{s)] a] ^ 
X]{n a„<a} f^n- As mentioned above this can be simply fixed by taking appropriate limits 
from the right. 

F<(a) = lim Cr\s-^P{s\A\. 

A trivial but useful comment is that F< is constant in intervals that do not contain points 
of A 

Example: A concrete example of the previous procedure is the following. Let us consider 
the sequences A = {an = n — 1 : n G N} and B = = 1 : n G N} associated with the 
generating function 

neN 

In this case 

F<{a)= I = 9{a){[a\ + 1) , a > 0. 

{neN:n-l<a} 

On the other hand 

neN 
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and it is clear that, for a non negative a ^ Nq, 



C-^ls-^P(sy.a]^-l_^^ (where ..>0) 



9{a - k) = 9{a){[a\ + 1) 

{feGNo:ifc<a} 



/ oo 
\ k=l 



sin 27rka 



where the last equality is obtained by using residues to compute the integral in the first line. 
As we can see in this case the inverse Laplace transform formula gives a + 1/2 for integer 
values of a (i.e. the average (F<(a + 0) + F<(a — 0))/2 of the left and right hmits). ■ 

Laplace transform and black hole entropy: The scheme presented above can be used 
compute 

iV<(a) = £-^[5-^ J] iV(a„)e-""^; a] , 

if a does not coincide with any of the values corresponding to the spectrum of the area 
operator, and extend the previous formula to all the positive values of the area -according 
to the entropy definition that we have adopted- just by taking limits from the right 

N^{an) = lim^ L'^ls'^ ^ A^(a,)e-"^^ a], a„ e Am- 

Notice that for N{an) > the previous limits are always well defined. 

The key point now is to realize that by using our generating functions, in particular 

G{l;xi,X2,---), we can get a simple expression for X^neN -^(^w)^"""*- ^^^^ '^^^ 
enough to substitute the arguments Xi in G{l;xi,X2,...) by Xi = e~*^/^. This is so 

because xf ■ ■ ■x'^'' ^ g-f (siv^+'-'+^r-v^) = q-o-s -v^^hen 2a = qi^/pi + h qr\fv'r- By doing 

this we find 

oo oo ^ 

P(s) := J2 ^^(«n)e"""' + 1 = ^(1; e-"v^/^ e-"v^/^ ■ ■ ■ ) = (l - 2 ^1 e~"^"^) ~ . 

neN i=l a=l 

The exponentials e~*^«^ appearing in this function can be simphfied by taking into account 
that the numbers are solutions to the Pell equation and hence 

This way ge get 

oo oo oo ^ 

i=l a=l k=l 

where we have used the fact that the values of the /c's appearing in the solutions to the 
Pell equations corresponding to different squarefree integers Pi are always different (so that 
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{/c^ : a G N} n {A;^ : a G N} = whenever i ^ j) and also that every /c G N appears in the 
solution to some Pell equation because can always be written as the product of 

a positive integer times a SRSFN. 

In order to take into account the projection constraint it is convenient to take z — e*'^/^, 
that in practice lets us get N{a,p) by performing an integral around a contour in the complex 
2;-planc consisting of a unit circumference surrounding the origin (notice that in this case 
z = 1 can be obtained by choosing w = 0). By doing this we get the function 



oo oo 



P{s, Lu) := G(e*^/2. _ _ _ {^~Y1 J^^e^''"^^ + e-^^^^/^^e-^^yJ^v^^ ^ . 



1=1 a=l 



The exponentials e "^y^VPi appearing in this function can be simplified if we use the Pell 
equations as before so we get 



oo oo 7 i 1 ^ 7 1 

P{s, uj)={l-2j2Yl e-\/^S(^/2 cos " = (l - 2 ^ e-\A(^/2 cos ^) " 

■ =1 a=l fc=l 



Notice that P(s, 0) = P{s). Finally, by performing the sums as explained above we have 
that 

s-'P{s,u) = s-^(l-2 5^e-^V^(^/2cos— ) 

k=i 

= P<{s,uj) + s-' (3.2) 

gives us the Laplace- Fourier transform P<(s,a;) of N<{a,p) plus a extra term that 
originates in the additional 1 appearing in the Domagala-Lewandowski prescription for 
n(a) = N<{a, 0) + 1. By inverting these expressions we then get the following result 

Theorem III.l. When a ^ Am, the values of N<{a), N<{a,p), and n{a) are given by 



N<{a) = / (s-^P(s,0) -s-Me»Ms 

2711 J^g_i^ V / 

= — / s-^(l-2j2 e-'V^)/A " ds - e(a) . 

1 /-Itt />x'()+ioo 

N<{a,p) = ^ / / (s-^P{s,u) - s-^Y'^'e-'^ dsdu 

= / / ,-i(i_2Ve-^v/^/^cos^:^) e^'e-'^dsdu 



-5(p,0)^(a). 

"Att pxo+ioo 



n{a) = / / s-^P{s,uj)e'''dsdu 

= / / s-i(i_2Ve-V^/2cos^) e-dsdu;, a > 0, 

oTT Z Jo Jxo-ioo ^ 2 / 
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where the value > is chosen to the right of all the singularities of the integrand in order 
to guarantee that the previous integrals converge. On the other hand, for a„ € Aih, the 
values A^<(a„), N<{an,p), and n(a„) coincide with the lim^_^^+ of the above expressions. □ 

At this point we should compare these results with the ones obtained by Meissner j^. 
First of all we see that once the extra 1 that appears in the definition of the entropy is 
incorporated the expression that we find for N<{a) exactly coincides with the one given by 
him. However, the expression that we find in the case when the projection constraint is 
taken into account differs from the one that appears in The difference between both 
expressions -once the extra 1 is taken into account- just amounts to an extra cosine factor 
in formula (35) of that should have been (using our P< for Meissner's P) 

P<is,u;) = 



2 °° / ujk\ / °° 

( exp{-s^k{k + 2)/4) cos j 1 - 2 ^ exp{-s^k{k + 2)/4) cos ^ 

^ k=l ^ ^ V fe=l 



Though this is somewhat speculative, there are two possible sources for this discrepancy: 
a simple typographical error or, more likely, an artifact introduced by the approximations 
that, as Meissner himself acknowledges, have been used in [3| to derive his formula (35). In 
our opinion the only way to settle this issue is to revisit Meissner's derivation and get his 
formulas again. This is the purpose of the next section. 



IV. FUNCTIONAL EQUATIONS FOR N<{a,p) AND N<{a) 

The main goal of this section is to obtain the results of the previous one by using functional 
equations as in [S]. We start by first considering the obtention of N^{a). The key idea is 
to pick a value for a and relate the values of N<{a) to those corresponding to a' < a. Let 
us consider the sequences in A/<(a) and classify them according to their first element. This 
allows us to partition this set disjoint union 

Ar<(a) = IjATf («) 

fceN 

with 

n 

Ar^^\a) := {m\3neN : me (Z,/2)" , ^ ^/\mi\{\mi\ + 1) < a, \mi\ = k/2} . 

Notice that Af<\a) = and N^^\a) = if a < ^3/2. The cardinalities of 7V<(a) and 
■^<\(^) (denoted by N^\a)) are then related by 

oo 

iV<(a) = 5^Arf (a). 

k=l 
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Notice also that the previous union involves, in practice, only a finite number of nonempty 
a) sets^° because given a fixed real value of a > the maximum of A; = 2|mi| for 

a sequence in Af<{a) is /cmax := [vl + 4a^ — Ij . 

For every value of /c < /cmax, and as long as a > there are two sequences in 

a) consisting of a single element, viz, {±k/2). The rest of them, if they exist, are of 

the form m = {±k/2, m2, . . .) and have, at least, two elements. If for each of the latter one 

considers the finite sequence = (777-2, • • •) the condition that m G equivalent 

to m, e^^<{a- ^k{k + 2)/2) so we conclude that the cardinality of each of the Af^^ (a) is 
simply given by 

N^^\a) = 2 + 2N<{a- ^k{k + 2)/2). (4.1) 

At this point it is convenient to extend the definition of N<{a) to arbitrary real values of a 
in such a way that N^{a) = for a < 0. This condition can be conveniently encoded as^^ 

N<{a) = e{a - V3/2)N<{a) , (4.2) 

and allows ( 14. ip to account for the case in which the sets N^\a) consist of sequences with 
a single element. 

Adding up the values of N^\a) given by equation (14.11) we get 

N<{a) = J2^<\^) = 2[Vl + 4a2-lJ+2^A'<(a- VA;(A; + 2)/2), a > ^3/2 . (4.3) 

k=l k=l 

Equation (14.31) for the function A^< : M No can be conveniently rewritten as 

00 

A^< (a) = 2 [V4a2 + 1 - Ij e{a - ^3/2) + 2 ^ A/'< (a - ^k{k + 2)/2) , (4.4) 

k=l 

where (14. 2 p allows us to extend the sum to infinity and the 9{a — a/3/2) factor is needed 
in the first term or the right hand side of (14. 4p to guarantee that it is zero for arbitrary 
negative values of a. Several comments are in order now. First of all it must be pointed out 
that formulas (14. 4p and (14. 2p correspond to equation (5) of Meissner ^ . We have carefully 
avoided to include the factor 6{a — V3/2) in (14. 4p as in [3| because it will be very useful 
to do so when deriving similar functional equations for N<{a,p). Notice also that condition 
(14. 2 p can be substituted by one of the type 

N<{a) =9{a- ao)N<{a) 

with < < VS/2 that is basically equivalent to requiring that N<(a) = for a < 0. 

We solve now the previous functional equations by using Laplace transforms. Although 
here we will just reproduce the correct equation (13) of Meissner, we give some details that 
will be relevant when discussing the resolution of the functional equations for N<{a,p). 
As already stated in ^ the fact that N<{a) is exponentially bounded j2| and piecewise 



For the same reason the previous sum is finite. 

Remember that with our definition for the function 9 we have 9{0) = 1. 
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continuous guarantees that its Laplace transform exists and is well defined in a half-plane 
{s G C : Re(s) > Xq} for some Xq G M. We have then 

P<(s) := [ N<{a)e-'''da (4.5) 

/» oo 
[V4a2 + 1 - 11 e-"' da + 2 f V A^< (a - v/fc(A; + 2)/2)) e""" da 



oo oo „ 

V- e-VMfc+2)/2 + 2 V / iV< (a - Vfc(fc + 2)/2) e"'^^ da 

oo oo „ 

^g-sy^^/2^2^e-^v/^(^/2 / iV<(a)e-'^Ma 

fc=l k=l 7[_^fc(fc+2)/2,oo) 

oo oo 

^ g-s^fc(fc+2)/2 _^ (s) ^ g-s^/fc(fcT2)/2 _ 



fc=l k=l 



Here we can change the order between sums and integrations in the two integrals appearing 
in the second line as a consequence of the Beppo Levi theorem^^. We have also used the 
fact that A^< (a) = for a < to set the lower limits in the integrals equal to zero in the last 
but one line of (14.51) . In fact this is the only condition that we need to solve the functional 
equation (^i^. We finally get jl] 



The fact that P<{s) is a proper Laplace transform tells us that we can write 

N<{a) = lim / P<{s)e^'ds. (4.7) 

2m A^a+ J^^_i^ 

for some a;o G M chosen in such a way that the singularities in the integrand are to the left of 
the integration contour in (14.71) . As we can see we recover precisely the same result obtained 
in the previous section by using our generating functions. 

A functional equation for N^{a,p) can be obtained in a similar way. Again we classify 
the sequences in J\f<{a,p) according to the first element and partition this set as the disjoint 



The Beppo Levi theorem is a corohary of the monotonous convergence theorem for Lebesgue integrals 
and states that if {/n}nGN is a sequence of non-negative measurable functions then 



oo p oo 
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union 

N<{a,p)= UAAf (a,p) 

with 



Af^\a,p) := {m|3neN : me (Z,/2)" , ^Tn^jf^n^TT) < a , = p , mi = A;/2} 



j=l i=l 



These are empty sets if a < \/3/2 and also if a < + 1). As before the cardinahty of 

J\f<{a,p) can be written in terms of N!^\a,p) as 



7V<(a,p)= J]7Vf (a,p), 



where the previous sum is, again, finite. 

Let us suppose now that p and a > \^/2. li k = 2p the only sequence of length 

one belonging to N'^\a,p) is (p); furthermore this can only happen if + 1) < a- 

If they exist, the remaining sequences in J\f!^\a,p) have, at least, two elements and are of 
the form m = (A';/2,m2, . . .)• We have now that m e J\f^\a,p) if and only if (m2, . . .) G 
Af<{a - ^k{k + 2)/2,p- k/2). li k ^ 2p there are no one-element sequences and the rest 
of them have, again, the form m = {k/2, 1712, ■ ■ ■)■ We can then write 

7Vf(a,p) = 5(fc,2p)%- vW^Tl)) + iV<(a- vW^T2)/2,p-fc/2), p^O. 

Here, as before, it is useful to extend the definition of A^< (a, p) to M and take this into account 
by imposing the condition A^< (a, p) = for a < 0. The previous reasoning is essentially valid 
in the p = case, the only difference is that, as the elements of the sequences arc non-zero 
half integers, it is impossible now to have unit length sequences. Summarizing we find that 

iVf (a,p) = (1 - 6i0,pmk, 2p)9{a - y/\p\{\p\ + 1)) + N<{a - V\k\m + 2)/2,p - k/2) . 
Adding up for all the possible values of k we get 

^<(«,P) = ^iVf (a,p) 



= (1 - 5(0,p))e(a - + + 5^ iV<(a - Vl^KI^I + 2)/2,p - k/2) 

= (1 - 5{O,p))0{a - V|p|(bl + 1)) (4.8) 

00 

+ J](A^<(a-y£(ZT2)/2,p-£/2) + iV<(a- v/^(7T2)/2,p + £/2)) . 
e=i 

Notice that it is immediate to check that by summing in p e Z/2 one recovers the functional 
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equation fl4.4p from (14.81) because 

5^ (l - 6{0,p))eia - VH(b| + l)) = 2[Va^^ - ma - V3/2) . 

It may seem a little bit surprising that we are not including another condition explicitly 
stating that N<{a,p) = if a < + 1) as done in jl]. In fact it is easy to see that this 

is a consequence of the functional relation (14 .Sp and the condition N<{a,p) = for a < 
just by repeatedly using it to compute N<{a,p) from values corresponding to smaller p's. 
The fact that we do not need to include a prefactor in the recurrence relation (14. 8p is the 
reason why it is indeed possible to get an exact solution to this functional equation without 
having to use any approximation. In our opinion the statement appearing in 13] claiming 
that one has to use approximations to solve the functional relation for A^< (a, p) stems from 
the difficulties in dealing with that prefactor. Let us solve (14. 8 p subject to the condition 
N^{a,p) = for a < 0. First notice that the fact that we can write 

N<{a)= N<{a,p) 
pez/2 

with N<{a) exponentially bounded implies that J2p(^z/2^^'^^-^<(^'P) exponentially 
bounded because 

and, hence, the previous sum has a well defined Laplace transform. The resolution of (14. 8 p 
is carried out as follows. 



[0,oo) 



J2 e'''^N^{a,p) 

P&Z/2 



da 



[0,oo) 



^^sal J2 e-^(l-5(0,p))%- v^HTHTI)) 



da 



pGZ/2 



[0,oo) 



pez/2 me 



da 



s ^ 2 

fc=i pez/2 



+5^6-^/ e-^'^^ iV<(a- v/|m|(|m|+2)/2,p-^)d 
pez/2 "^[0'°°) mez. 

oo , 

^ g-Vfc(fc+2)/2 cos + ^ gi^P ^ g-s^|m|(|m|+2)/2 



fc=l 

oo 



pez/2 

oo 



W<(a,p 



)da 



[0,oo) 



^ ^~s^k^)l2 cos + 2 5^ ^-s^Uk^)l2 ^ g.c.p' / g-.a^^ p/)^^ 

fc=l fc=l p'eZ/2 >^[0,oo) 

oo , oo , 

^ g-s V^(^/2 fc^^ ^ 2P< (s, o;) 5^ e-^v/^(^/2 fc^^^ _ 



fc=i fc=i 
Here we can justify to change the order of sums and integrals as a consequence of the 
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dominated convergence theorem. We have also used the fact that A^< (a, p) = for a < in 
order to write the lower limits in the previous integrals equal to zero as we did for P<{s). It 
is important to notice at this point that this is the only condition that we need to impose in 
order to solve the previous functional relations. In particular we have been able to do this 
without the 6{a — a/|p|(|p| + 1) prefactor used by Meissner. In our opinion this is crucial to 
avoid the use of simplifying assumptions. According to the previous derivation we get 

P< {s, 00) = -^(^^ e~^V^)/^ cos ^ j ^1 - 2 e-^^^y^ cos ^ j , (4.9) 

in perfect agreement with the result obtained in the previous sections by using generating 
functions. We can get N<{a,p) from (14. 9 p by inverting the Laplace- Fourier transform and 
recover the result of theorem llll.li 



V. ANALYTIC PROPERTIES OF P{s) 

We study in this section the analytic structure of the function P{s). This is a first 
necessary step to understand the behaviour of P{s, u). We start by enunciating and proving 
several lemmas concerning the poles of P{s) and finally give the main result of this section 
concerning the accumulation of the real parts of these poles precisely to the value jm- This 
result is important in the asymptotic analysis of the black hole entropy. In the following we 
will denote 

^ oo 

Q(s) :=-±- = l-2j2 e-'V^^)l^ , Re(s) > . 
^^■'^^ k=i 

This is obviously an analytic function in any band Re(s) > a;o > due to the uniform 
convergence of the series 

Er=i e-Vfc('=+2)/2. In fact 



Vfc(fc+2)/2| ^ g-Re(s)^fc(fc+2)/2 ^ ^-xo^ k{k+2) /2 



for all s in the strip Re(s) > xq > 0. Hence, using Weierstrass criterion, Yl^=i e~*v'*^^*^+2)/2 
converges uniformly in Re(s) > xq > to an analytic function (see, for example, [13] 
for a good review of convergence properties of Dirichlet series). In fact, it can be proved 
that Q is an analytic almost periodic function [l3|. Remember that an analytic function 
/(s) = /(x + iy), regular in a strip xi < x < X2 (— oo < Xi < X2 < +oo), is called almost 
periodic (uniformly almost periodic) if for e very e there exists a length L = L{e) such that 
every interval yo < y < i/q + L of length L on the imaginary axis contains at least one 
translational number r = T{e) associated with number r satisfying the inequality 

\f{s + zT)-f{s)\<e 



for all s in the strip xi < x < X2. Every periodic function, such as fk{s) = e ^v^(^+2)/2^ jg 
almost periodic. Also a uniformly convergent sequence of almost periodic functions, such as 
Er=i e-V^('=+2)/2^ is almost periodic [13. 



Prior to giving the main result of this section we state and prove some simple but useful 
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lemmas. 

Lemma V.l. Let F be the restriction of Q to the positive real axis = (0, oo), then F is 
an analytic and monotonically growing function. 

This is so because 1 — F is the hmit of a sum of strictly monotonically decreasing functions 
of the type e-^V^^^+^vi. □ 

Lemma V.2. F has a single zero in and, hence, Q has only one real zero. 

This is an immediate consequence of the continuity and monotonicity of F and the fact that 



lim F{x) = — oo , lim F{x) = 1 . 



2 

k=l 



Lemma V.3. The only zero of Q with real part equal to 'Jm is 7a/ itself. 
Let us suppose that there exists sq = 7m + iyo G C, with yo such that 



fc=l 

Then the following conditions must hold 



k=l 

oo 



k=l 



However, from ( I5.ip we must have 



^e-7My^(F^/2 ^ J2 e-'^M^m^)l2 cos{yo,/k{k + 2)/2) 



k=l k=l 

which is impossible. In fact 



□ 



In the following we will denote this zero as 7a/. It obviously satisfies 

oo ^ 



oo ^ 

Yl e-^^'V^)/^ cos(^oVMfc + 2)/2) = - , (5.2) 



e-7MVfc(fc+2)/2 sm{yo^/k{k + 2)/2) = . (5.3) 



e-7MVMfc+2)/2 > g-7MVMfc+2)/2 cos{yo^k{k + 2)/2) , VA; G N 
and taking into account the fact that ^/k{k + 2) is irrational^^ and yo ^ 0, there exists 



If yo\Jk{k + 2)/2 = mod 27r for some k we can always find kg ^ Z such that yoy/ko{kQ + 2)/2 ^ 
mod 27r. 
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fcn G N such that 



e-7A/V*^o(fco+2)/2 > g-7M Vfco(fco+2)/2 cos(yo ^^^0(^:0 + 2) /2) 

We then conclude that 



^e-7M 7^(^/2 > ^ g-7M 7^(^/2 cos(yo VA;(A; + 2) /2) , 
fc=i fe=i 

and, hence, (15. 2p cannot be satisfied. We then conclude that sq cannot be a zero of Q. □ 

Lemma V.4. The real part of the zeroes of Q different from 'is strictly smaller than ^m- 
Let us take Sq = Xq + iijo G C, Sq 7^ 7m, such that Q{so) = 0. Then it follows that 



1 
2 



k=l 



00 
k=l 



^g-xoyfc(fc+2)/2 _ 



Using now (15.11) . we find that 

00 00 

^ ^~jMy/k{k+2)/2 ^ ^ ^~Xoy/k{k+2)/2 _ ^^^^^ 



k=l k=l 



Now, taking into account that the function / : (0, 00) ^ M : x J2kLi e-'^v^^^/^ is 
continuous and monotonically decreasing, we see that equation (15.41) implies that 

As we have assumed that sq 7^ we conclude that Re(so) = xq < 7m- D 
We end this section with the following important result 

Theorem V.5. The function Q has an infinite number of zeros. The set of the real parts 
of the zeros of Q has an accumulation point at 7m- 

An immediate consequence of almost periodicity and uniform convergence in a strip is 
that, if the equation Q{s) = is solvable in the strip Re(s) > Xq > 0, then it will have 
infinitely many solutions and their imaginary parts will form a relatively dense set^"^. This 
is a well-known application of Rouche's Theorem. Here we basically repeat the reasoning 
appearing in [l3|. Let 7m be the real zero of Q{s) = in the strip Re(s) > xq > 0. Then 
there exists an ro > such that the circumference C;y^j-(ro) = {s G C : |s — 7m | = '''o} is 



A set T of real numbers r is called relatively dense [l3| if there are no arbitrary gaps among the numbers 
r or, to be exact, if some length L exists such that every interval {a,a + L) C M of this length contains at 
least one number r of the set T. Roughly speaking, a relatively dense set can be described as one that is 
just as dense as an arithmetic progression {an : n e Z}, a > 0. 
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contained in the strip, encircles a single zero of Q{s), and Q{s) 7^ on the points of C;y^^(ro). 
Take now 

e := min{|g(s)| : s G (^^^(^o)} • 

By uniform almost periodicity, there exists L > such that every interval of length L 
contains r G M such that \Q{s + it) — Q{s)\ < e along C;y^^(ro) and, hence, 

\Q{s + it) - Q{s)\ < \Q{s)\ , Vs G C^,,(ro) . 

By Rouche's Theorem, we deduce that Q{s) and Q{s + it) have the same number of zeros 
inside the circle C^^^^(ro). By repeating the argument for every < r < tq and taking into 
account that all zeros of Q (except 7j\/) satisfy Re(s) < 7m, we see that we can find zeros 
of Q with real parts smaller but as close to 7Af as we wish. □ 

As shown in the following example (which is a simple extension of the example given 
above and of the same type as those given by Meissner in [sl) the accumulation of the real 
parts of the poles in the integrand of an inverse Laplace transform changes its asymptotic 
behavior relative to the one that one would expect by considering only the real pole^^. 

Example: Let us consider now the sequences A = {an = n — 1 : n G N} and B = = 
2n-i ■ 7^ g pjj associated with the generating function 

G{x) = y{2x''-' ^ 



- 1 - 2x • 

neN 

In this case 

F<(a) = Yl 2"-i = ^(a)(2L'^J+^-l) 

{nGN:n-l<a} 

On the other hand 



-(n-l)s 



1 - 2e~ 



and we have now for a non negative a ^ Nq, 

I i-XQ+ioo ^as 

C-^[s-^P(sy,a] = — (where Xo>0) 

^T^^Jxo-ioo s{l-2e 

J2 2^=0(0 - /t) = ^(a)(2L"J+i-l) 

{fceNo :fc<a} 

AkiT sm{2ak7i) + (2 log 2) cos(2aA;7r) 



log 2 ^ 47r2A;2 + log^ 2 

where the last equality can be obtained by using residues to compute the integral. We can 
see in this case that the inverse Laplace transform formula gives 3 ■ 2*^*^^^^ — 1 for integer 



A small perturbation of the previous function would allow us to have accumulating real parts different 
from the value corresponding to the real pole in the integrand. 
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values of a (as expected the average (F<(a + 0) + F<(a — 0))/2 of the left and right limits). 
The important issue now is to realize that 2" — 2^*^^ oscillates with an exponentially growing 
amplitude and hence the values of F<{a) for a ^ oo are not simply proportional to 2". ■ 



VI. CONCLUSIONS 

In this paper we have discussed several ways to exactly compute the black hole entropy 
in loop quantum gravity [ij. In particular we have given a procedure based on generating 
functions that gives an independent way to derive the results appearing in the literature on 
this issue. In this way we have been able to detect and correct a mistake in the expression 
appearing in [st] for the black hole entropy when the projection constraint is taken into 
account. Second we have shown that the number-theoretical methods introduced in 0, 0], 
and successfully used to get precise numerical information about the entropy for small black 
holes, can be used in an efficient way to obtain exact formulas in the spirit of f^. In a sense, 
the two approaches are unified in this paper. We would like to point out, anyway, that 
the kind of detailed information provided by the number-theoretical methods of [3] is very 
difficult to extract form the expressions of the entropy as integral transforms. Finally, we 
have discussed the analytic properties of some of the functions that appear in the expressions 
of the black hole entropy. The most important result in this respect concerns the distribution 
of poles in the integrand of the inverse Laplace transform defining the entropy. From this 
analysis we have shown that the value of the Immirzi parameter given in the literature is, 
in a sense, correct but the asymptotic analysis of the entropy may display an interesting 
behaviour superimposed on the expected linear growth of the entropy as a function of the 
area. This suggests that the entropy structure found in 0] for small black holes in 
numerical computations may actually survive for macroscopic areas. We want to add several 
comments: 

In our opinion the paper by Meissner is usually misread and misunderstood. In particular, 
the exponential form of the entropy in terms of the area is often taken as some kind of 
ansatz introduced to approximately solve the functional equations giving the black hole 
entropy (with and without the projection constraint). It is very important to emphasize 
that by using Laplace transforms or Fourier-Laplace transforms it is possible to get exact 
expressions for the entropy. The exponential ansatz can be used to quickly show that the 
sought exponential growth of the entropy is somehow present but should not be taken as a 
rigorous derivation of the asymptotic behaviour of the entropy. 

Despite the claims by the author of it is possible to find an exact closed expression for 
the black hole entropy also when the projection constraint is included. We have shown this 
in two complementary ways. By using the generating functions given in Q and by solving a 
functional equation along the lines suggested in [3]. In the latter case we have corrected an 
error in the original functional equation and rewritten it in a way that facilitates its exact 
resolution. Actually the solution given in is very close to the right one. 

We have studied in detail the analytic structure of the functions appearing in the in- 
tegrand of the contour integral that gives the entropy (in the simplified setting where the 
projection constraint is not incorporated) to see if its asymptotic expansion can be readily 
obtained by looking at the poles of the integrand as claimed in ^ . We have proved several 
important results in this respect (in some cases completing the claims of Meissner): 



i) There is indeed an infinite number of poles in the integrand of (14. 5p . 
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ii) They are confined to a band in the complex plane and their real parts are bounded 
from above by ^m- 

iii) There is only a single pole of the integrand with real part equal to 7m- 

iv) The real parts of the poles have an accumulation point precisely for the value 7m (and 
maybe others). 

The last point is specially important as far as the asymptotic behaviour of the entropy is 
concerned because in a situation such as the one described here it may not be true that the 
asymptotic behaviour is given by the contributions of the pole with the largest real part, in 
fact this is illustrated by the simple example given in section |Vl Notice also that even if this 
is the case the fact that one has other poles with real parts arbitrarily close to 7m means 
that there may be corrections that are relevant for macroscopic but not infinite areas. It is 
clear, nonetheless, that the exponential behaviour given in [3j, should play a significant role 
in the final asymptotic form of the entropy. In particular the value of the Immirzi parameter 
7Af is distinguished by the fact that is the maximum of the real parts of the poles. In this 
sense it controls the growth of the entropy despite the fact that its asymptotic behaviour 
may not be given by a simple exponential. 

The reader may argue that, in the end, the methods based on the solution of functional 
equations for the entropy as introduced by Meissner are rather quick and efficient so there is 
no need to resort to the kind of detailed combinatorial analysis of 0, 0]. Our opinion is that, 
although they are indeed very clever ways of tackling this hard problem, the kind of detailed 
information provided by our number-theoretical approach is very useful as shown by the fact 
that they provide an independent way to check the results obtained so far. Taken at face 
value the expressions for the black hole entropy as inverse integral transforms give the entropy 
at all scales (and hence also in the microscopic regime where the interesting behaviour of the 
entropy has been found [^). However, due to the subtly oscillatory nature of the integrands 
they cannot be practically used to obtain the entropy with good precision. Our combinatorial 
methods are much better in this respect and, in any case it is possible to check that they 
give the same values for the entropy that the integral expressions (whenever they can be 
numerically computed). Also, they allow to exactly characterize the area spectrum and the 
microscopic configurations corresponding to the allowed values of the area going far beyond 
the results obtained in 

The issue of getting the right asymptotic and the behaviour of the entropy for macroscopic 
scales, for which the integral expressions that we give here are a good starting point, is the 
last important problem that remains to be addressed to completely understand the behaviour 
of the black hole entropy in LQG. This will be our goal in the immediate future. 
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